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Abstract 

Applying the perturbative QCD ( PQCD ), we study the process B c — > D*^ in 
the technicolor with a massless scalar doublet model ( TCMLSM ) recently presented 
by C. D. Carone and H. Georgi; and compare the results with that estimated in 
the standard model ( SM ). There are two mechanisms which contribute to the 
B c — ► D*j process. One proceeds through the short distance b — ► S7 transition, the 
other through weak annihilation accompanied by a photon emission. In the SM, 
these two processes are found to contribute with the same order of magnitude. In 
the TCMLSM, the modification of B c — * D*^y from tt p ( the physical pions in this 
model ) is small for the allowed mass of tt p . 

PACS numbers: 12.15.LK, 12.60.Nz, 13.30.Eg 
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I. Introduction 



The inclusive rare decay B — > X s ^ has been studied several years before [1] . Recently 
the physics of B c meson has caught intensive attentions [2]. The B c meson is believed to 
be the next and the final family of B mesons. It provides unique opportunity to examine 
various heavy quark fragmentation models, heavy quark spin-flavor symmetry, different 
quarkonium bound state models and properties of inclusive decay channels. Being made 
of two heavy quarks of different flavors, B c radiative weak decay also offer a rich source 
to measure elements of Cabibbo-Kobayashi- Maskawa ( CKM ) matrix of the standard 
model ( SM ). 

Different from the general rare B decays B — > X s 7 which is mainly induced by the 
flavor-changing b — * 57 neutral currents [3], in B c — > D*7, the bound state effects could 
seriously modify the results from the assumption. Bound state effects include modifi- 
cations from weak annihilation which involve no neutral flavor- changing currents at all. 
The effects of weak annihilation mechanism are expected large due to the large CKM 
amplitude. 

Unfortunately, the well known chiral-symmetry [4] and the heavy quark symmetry 
[5] can not be applied to B c — > D*^y process. Recently, a perturbative QCD ( PQCD 
) analysis of B meson decays seems to give a good prediction [6]. As it is argued in 
Ref . [7] , B c two body nonleptonic decay can be conveniently studied within the framework 
of PQCD suggested by Brodsky-Lepage [8] and then developed in Ref. [6]. Here, we 
summarize their idea: In the subprocess b — > S7, s quark obtains a large momentum by 
recoiling, in order to form a bound state with the spectator c quark, the most momentum 
of s must be transferred to c by a hard scattering process. In the final bound state ( i.e. 
D* ), since the heavy charm should share the most momentum of D* s , the hard scattering 
is suitable for PQCD calculation [6, 8]. 

In B c — > -D*7, the subprocess b — > sj, taken as a free decay, is usually controlled 
by the one-loop electromagnetic penguin diagrams which are in particular sensitive to 
contributions from those new physics beyond the SM. This situation is similar with B — > 

Recently, the modification of the electromagnetic penguin interaction to the decay 
b — > 57 from PGBs in the one generation technicolor model ( OGTM ) has been estimated 
in Ref. [9]. In the recent literatures [10, 11], the technicolor ( TC ) models with scalars 

2 



have been studied extensively. The phenomenological studies have shown that the TC 
models with scalars do not generate unacceptably large FCNCs and are consistent with 
the experimental constraints on oblique electroweak radiative corrections. Among these 
models, the TC with a massless scalar doublet model ( TCMLSM ) [10], presented by C. 
D. Carone and H. Georgi is the simplest nontrivial extension of the SM with only two 
new free parameters ( h + , h_ or h, A ). The phenomenology of the TCMLSM has been 
discussed in the literatures [11]. In this paper we express these two new free parameters as 
m np and jj differently from Ref. [10]. The relation of them can be found in Ref.[10], where 
m np is the mass of i\ v ( the physical pions in this model ), / and /' are the technipion 
decay constant and scalar vacuum expectation value ( VEV ), respectively. The couplings 
of charged n p with ordinary fermions are given as 

\/2 f 

K - -Ui - -dj] = — —V Uidj — [m dj (l + 75) - m Ui (l - 75 )], (1) 

y/2 f 

K - ~ u i - -dj) = — ^-K^-^KJl + 75) - m dj (l - 75 )], (2) 

where u=( u, c, t ), d=( d, s, b ), V Uidj is the element of CKM matrix, and v ~ 250 GeV 
is the electroweak scale. 

In this paper, applying the above PQCD method, we address B c — > _D*7 in the 
TCMLSM to examine the virtual effects of n p in the TCMLSM and compare the results 
with which estimated in the SM. 

This paper is organized as follows: In Sec. II, we display our calculations in the SM 
and TCMLSM. We present the final numerical results in Sec. III. Sec. IV contains the 
discussion. 



II. Calculation 

Using the factorization scheme [8] within PQCD, the momentum of quarks are taken 
as some fractions x of the total momentum of the meson weighted by a soft physics 
distribution functions $ H (x). The meson wave functions of B c and D* s take the simple 
form of 5 function ( the so-called peaking approximation ) [12, 13]: 

= 2v /3 > ( 3 ) 
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= ,D; % eD; \ (4) 



The normalization [13] is 



Jo 



dx® Bc {x) = (5) 







where x, y denote the momentum fractions of c, s quarks in the B c and Z)* mesons, fs c 
and /dj are decay constants of B c and D* respectively, 

= m D .-m c 

The spinor part of B c and D* [14] are 

(j? + ™B c )75 



V2 ' 



(9) 



(10) 



V2 ' 

which come from the matrix structures of B c and D* meson wave functions, while p and 
q are the momenta of the B c and D* respectively, and e is the polarization vector of D* s . 

II ( i ). Electromagnetic penguin contribution 

The relevant Feynman diagrams which contribute to the short distance electromagnetic 
penguin process b — > 57 are illustrated as the blob of Fig.l. 

In the evaluation, we at first integrate out the top quark and the weak W bosons at 
H = mw scale, generating an effective five-quark theory. By using the renormalization 
group equation, we run the effective field theory down to b-quark scale to give the leading 
log QCD corrections. 

After applying the full QCD equations of motion [15], a complete set of dimension-6 
operators relevant for b — > 57 decay can be chosen as 0\ - 0$, which have been given in 
the Refs.fl, 9]. The effective Hamiltonian appears at the W scale is given as 



AG 



Heff = -j=-V tb VZ'£C i (m w )O i (m w ). (11) 
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The coefficients of 8 operators are: 



Ci(m w ) =0,1 = 1,3, 4, 5, 6, C 2 (m w ) = -1, (12) 
C 7 (mw) = \A{x) - (jj) 2 [B(y) - l -A{y)], (13) 

Cs(m w ) = l -D{x) + {^) 2 \\D{y) - E(y)}, (14) 

where functions A, B, D and E are defined in the Refill, x = ( JB M 2 , y = f-^) 2 . 

The running of the coefficients of operators from \x = mw to /x = m b was well described 
in Ref. [16]. After renormalization group running we have the QCD corrected coefficients 
of operators at /i — m b scale. 

16 8 i 16 8 

C e 7 }f {m b ) = Q-%C 7 (m w ) + -(g-H - g -£)C 8 (m w ) + C 2 (m w ) ]T h t Q~ a % (15) 

d i=i 

with 

£ = ^ v (16) 

a s (m w ) 

,626126 56281 3 1 . . 

hi = (zmfj, -sTm* ~r "i4' ~ - 6494, _0 - 0380, "°- 0186, -°- 0057 )' 

ai = ( 2^' 2^' 4' °' 4 ° 86 ' -°' 4230 ' -°- 8994 ' °- 1456 )- ( 18 ) 
Now we write down the amplitude of Fig.l as 

M a = ^dx 1 dy 1 ^(y 1 )^x 1 )^V tb V*C e 7 ff (m b )m b e^C F 

{T r [(/i - m* D J jta^l + l5 )k"rf(j>- Vl A + m b ) la (i> + m Bc ) lbl ^_^_ (19) 

where rj is the polarization vector of photon, x\, y\ are the momentum fractions shared 
by charms in B c and D* respectively, and the color factor Cf — 3- The factors D±, D 2 
and -D3 in eq.(19) are the forms of 

Di = (1 -yi)(m 2 Bc -m 2 D , yi ) -m 2 b , (20) 
D 2 = (l-x 1 )(m 2 D .-m 2 B x 1 ), (21) 
D 3 = (xi - yi)(xim 2 Bc - yim 2 D *). (22) 

Now the amplitude M a can be written as 

M a = ie^k'eYfr 9 + V^Mk -m 2 D J-(p + q),(e ■ k)]fr 9 , (23) 
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with form factors 

fr 9 = Vl en9 = C f} dx 1 dy 1 6(x 1 - e Bc )5( yi - e D *) 

{[m Bc {l - Vi)(m Bc ~ 2m D ») - m b {2m Bc - m DJ )]^^ - m Bc m D *(l - x^^-J, 

' (24) 

where 

c = em b f Bc f D * C? ff (m b )C F a s (m b )G F V tb V t * 

12-\/2 

II ( ii ). The weak annihilation contribution 

As mentioned in Sec. I, B c meson is also the unique probe of the weak annihilation 
mechanism. The leading log QCD-corrected effective weak Hamiltonian of W annihilation 
is 

H$ = ^V cb V; s (c + + + c_0_) + H.c, (26) 

with 0± = (sb)(cc) ± (sc)(cb), where (q^) = <?i7/i(l — 75)^2, c± are Wilson coefficient 
functions. 

Using the method developed by H. Y. Cheng et al. [17], we can get the amplitude of 
W annihilation diagrams ( see Fig. 2 ): 

M b iW) = ie„ ua ^k»e a pPf?™ + rK(mk ~ m 2 Ds ) " (p + 9)„(c • WfSSfi, (27) 



with 



where 



Aw = m- + + (- + -)] T- mB i . ( 2 8) 

; m s m c m Bc m c m b m Bc — m D * 

/So = -CK^ - + - (29) 

^ ; m s m c m Bc m c m b m Bc — m D * 

( = ea 2 ^=V cb V;j B J Dh (30) 

and 02 = |(c_ — c+) is a calculable coefficient in the nonleptonic B decays. 

Using the Feynman rules given in eq. ( 1 ) and eq. ( 2 ), the leading log QCD-corrected 
effective weak Hamiltonian of annihilation is given as 

H { : f j = - V ^s^f^^ + + + c_0_) + H.c. (31) 
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Using the same method as the above, we can write down the amplitude of ir^ annihi- 
lation diagrams ( see Fig. 2 ) as 

M b M = * W r^eVA?C) + -"&;)- (P + ?),.( e ' k)]f%%, (32) 



with 



where 



/i"","") = 2C'[(- + ^)^^£ + (iS. + f^I^^^S., (33) 
k p; m s m c m Bc m b m c m Bc m% c - m z D , 

/S*) = C K- + + (- + * )]-?3l!!*_, (34) 

C' = e^K.KU^) 2 ^^/^^ + my- (35) 

*/ 7Tt> 



(36) 



The total annihilation amplitude ( Fig. 2 ) is the form of 

M b = M b {w) + M b M 

= ^a/^A^V/r™ + vn^« c -m 2 Ds )-{p + g) M (e • fc)]/ 2 on ™, 

where 

ranni ranni , ranni fQ f 7\ 

Jl — Jl(W) "+" A(7r p )) l" 3 ' J 

ranni ranni , ranni fQQ\ 

J 2 — J2(W) + /2(7rp)- I' 33 J 

Finally, we estimate another possible long-distance effect, namely the vector -meson- 
dominance ( VMD ) contribution which was advocated by Golowich and Pakvasa [18]. 
VMD implies that a possible contribution to B c — > D*^ comes from the B c — > D*J/ip(ip') 
followed by J/tp(ip') — * 7 conversion. 

As discussed in Refs.[19, 20], using factorization approach, the VMD amplitude ( Fig. 3 

)is 

M VMD = ie^TfkTfpPf™ + <[e„(m| c -m 2 Ds )-(p + q),{e ■ k)}g MD , (39) 

with 



fVMD =eG F V db V*{y/2a 2 1 ( jiZW^ + tofei) 

+a 1 f D .m D .[-, + 7 r 1 -^ ]}^ BD (0), 



(40) 



rVMD _ _ 1 D (Q) eVMD (a 1 \ 

^2 - n T /BD/n\Jl ' ^ 4i J 



2\/ BD (0) 

where ai = |(c+ + c_), and l /BD (0) and Af D (0) are form factors. 



7 



III. Numerical results 



We will use the following values for various quantities to carry on our calculations. 
( i ). Decay constants for mesons. Here we use 

fD- = fn s = 344MeV[21], f Bc = 500MeV[22], f J/f = 395MeV[20], fa = 293MeV[20]. 

( ii ). Meson mass and the constituent quark mass [23, 24] 

m Bc = 6.27GeV, m D . = 2.11GeV, m b = 4.7GeV, m c = 1.6GeV, 

m s = OMGeV, rrij/^ = 3.079GeV, = 3.685GeV. 

We also use m Bc m h + m c , m D * m s + m c in our calculations. 

( iii ). ai and a 2 have been estimated very recently in Ref.[17] according to the CLEO 
data [25] on B — > D*n(p) and B — > J/ipK*. Here we take 

ai = 1.01, 02 = 0.21. 

( iv ). CKM matrix elements [24]. Here we use 

|F cb | = 0.04, \V ts \ = \V cb \, \V CS \ = 0.9745, |V^| = 0.9991. 

( v ). The QCD coupling constant a s (/x) at any renormalization scale, can be calculated 
from a s {m z ) = 0.117 via 

= ^ 1 • 

1 -(11- 2^)2^1 M(^) 

We obtain 

a s {m b ) = 0.203, a s (m w ) = 0.119. 

( vi ). The Ref.[10] gives a constraint on m^ p in the allowed parameter space of the 
model: \m z < m Wp < ITeV. Here we take 

m np = (50 ~ 1000)Gey. 

( vii ). With the constraint of f 2 + f' 2 = v 2 and the chiral perturbation theory in 
Ref.[10], we can get 0.115 <j< 1-74. Here we take 

T = °' 115 
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in our calculations. 

( viii ). The form factors V(0) and A 2 (0) appearing in the two-body decays of B. 
From Ref. [26], we take 

V BD (0) = 0.30, Af D (0) = 0.20. 

We present the form factors ft ( /f e " 3 , /f ng , f? nn \ f£ nni ) in the SM and TCMLSM 
in Table 1, so do the decay widths in Table 2 using the amplitude formula 

The calculated results indicate that the VMD effects are large near the pole and which 
can not be neglected. The calculated results are 

f VMD = g 73 x 10 -1 0> f VMD = x 10 -1 0> 

t vmd = 112x w-^GeV. 
The lifetime of B c is given in Ref. [27]. In this paper we use 

t b c = (OAps ~ 1.35ps) 

to estimate the branching ratio BR ( B c — > _D*7 ) which is a function of t Bc . The results 
are given in Table 3. 



IV. Discussion 

Applying PQCD, we have studied two mechanisms which contribute to the process 
B c — > .0*7. For the short-distance one ( Fig.l ) induced by electromagnetic penguin 
diagrams, the momentum square of the hard scattering exchanged by gluon is about 
3.QGeV 2 which is large enough for PQCD analyzing. The hard scattering process can not 
be included conveniently in the soft hadronic process described by the wave-function of 
the final bound state, which is one important reason that we can not apply the commonly 
used models with spectator [28] to the two body B c decays. There is no phase-space for 
the propagators appearing in Fig. 1 to go on-shell, consequently, unlike the situation in the 
Ref. [6], the imaginary part of M a is absent. Another competitive mechanism is the weak 
annihilation. In the SM, we find this mechanism is as important as the former one ( they 
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can contribute with the same order of magnitude ). This situation is distinct from that 
of the radiative weak B ± decays which is overwhelmingly dominated by electromagnetic 
penguin. This is due to two reasons: one is that the compact size of B c meson enhances 
the importance of annihilation decays; the other comes from the Cabibbo allowance: in 
B c — > -D*7, the CKM amplitude of weak annihilation is IV^V^J, but in B ± — > K ± ^ ) the 
CKM part is |KbK*J which is much smaller than 

Particularly, the VMD contribution is found not small. This situation is quite different 
from the cases B — > J/ipK(K*) and B — > J/ipp [19, 20]. The reason comes from that 
although the coupling of J/ip(ip') — 7 is small ( e/g-yj/^') ~ 0.025(0.016) ), the J/ip(ip') 
resonance effect can be very large. 

In addition, we find that the modification of B c — > D*^ from ir p in the TCMLSM is 
small for the allowed range of mass of n p ( with -jj fixed ). This situation is quite different 
from that of Ref . [9] , in which the size of contribution to the rare decay of b — > S7 from the 
PGBs strongly depends on the values of the masses of the charged PGBs. The difference 
is mainly due to the small value of -jj which leads to the small modification from ir p in 
the TCMLSM. However, in the OGTM, such suppression factor jj does not exist. In 
our calculations, we take -p=0.115 as the input parameter. When -p is taken properly 
larger ( without exceeding the constraint: 0.115 < -jj < 1.74 ), the calculated results 
remain unchanged basically. In view of the above situation, it seems to indicate that the 
window of process B c — > D*^ is close for the TCMLSM. But in our calculations, besides 
the peaking approximation of the meson wave functions, the theoretical uncertainties 
are neglected, such as that of a(m z ), next-to-leading log QCD contribution [29], QCD 
correction from m t to m w [30], etc. When the more reliable estimation is available within 
the next few years, one can, in principle, make the final decision whether the window for 
TCMLSM is open or close. 

ACKNOWLEDGMENT 

This work was supported in part by the National Natural Science Foundation of China, 
and by the funds from Henan Science and Technology Committee. 

Reference 



10 



1. B. Grinstein et al, Nucl. Phys. B 339 ( 1990 ) 269. 

2. D. S. Du and Z. Wang, Phys. Rev. D 39 ( 1989 ) 1342; K. Cheng, T. C. Yuan, 
Phys. Lett. B 325 ( 1994 ) 481, Phys. Rev. D 48 ( 1994 ) 5049; G.R. Lu, Y.D. 
Yang and H.B. Li, Phys. Lett. B341 (1995)391, Phys. Rev. D51 (1995)2201. 

3. J. Tang, J. H. Liu and K. T. Chao, Phys. Rev. D51 (1995)3501; K. C. Bowler 
et al, Phys. Rev. Lett. 72 ( 1994 ) 1398. 

4. H. Leutwyler and M. Roos, Z. Phys. C 25 ( 1984 ) 91. 

5. M. Neubert, Phys. Rep. 245 ( 1994 ) 1398. 

6. A. Szczepaiak et al, Phys. Lett. B 243 ( 1990 ) 287; C. E. Carlson and J. Milana, 
Phys. Lett. B 301 ( 1993 ) 237, Phys. Rev. D 49 ( 1994 ) 5908, ibid 51 ( 95 ) 450; 
H-n. Li and H. L. Yu, Phys. Rev. Lett. 74 ( 1995 ) 4388. 



7. D. S. Du, et al, BIHEP-TH-95-38 ( |hep-ph/ 96032911 



8. S. J. Brodsky and G. P. Lepage, Phys. Rev. D 22 ( 1980 ) 2157. 

9. Cai-Dian Lii and Zhenjun Xiao, Phys. Rev. D 53 ( 1996 ) 2529. 

10. C. D. Carone and H. Georgi, Phys. Rev. D 49 ( 1994 ) 1427. 

11. C. D. Carone and E. H. Simmons, Nucl. Phys. B 397 ( 1993 ) 591; E. H. Simmons, 
Nucl. Phys. B 312 ( 1989 ) 253; S. Samuel, ibid, B 347 ( 1990 ) 625; A. Kagan and 
S. Samuel, Phys. Lett. B 252 ( 1990 ) 605, B 270 ( 1991 ) 37. 

12. S. J. Brodsky and C. R. Ji, Phys. Rev. Lett. 55 ( 1985 ) 2257. 

13. Hsiang-nan Li, Phys. Rev. D 52, ( 1995 ) 3958. 

14. C. E. Carlson, J. Milana, Phys. Rev. D 51, ( 1995 ) 4950. 

15. H. D. Politzer, Nucl. Phys. B 172 ( 1980 ) 349; H. Simma, Preprint, DESY 93-083. 

16. M. Misiak, Phys. Lett. B 269 ( 1991 ) 161; K. Adel, Y. P. Yao, Mod. Phys. Lett. 
A 8 ( 1993 ) 1679; Phys. Rev. D 49 ( 1994 ) 4945; M. Ciuchini et al, Phys. Lett. 
B 316 ( 1993 ) 127. 



11 



17. H. Y. Cheng et al., Phys. Rev. D 51 ( 1995 ) 1199. 

18. E. Golowich and S. Pakvasa, Phys. Lett. B 205 ( 1988 ) 393. 

19. M. Bauer, B. Stech and M. Wirbel, Z. Phys.C 34 ( 1987 ) 103; E. Golowich and S. 
Pakvasa, Phys. Rev. D 51 ( 1995 ) 1215. 

20. H. Y. Cheng, Phys. Rev. D 51 ( 1995 ) 6228. 

21. A. Aoki et al., Prog. Theor. Phys. 89 ( 1993 ) 137; D. Acosta et al, CLNS 
93/1238; J. Z. Bai et al., BES Collaboration Phys. Rev. Lett. 74 ( 1995 ) 4599. 

22. W. Buchmiiiller and S-H. HTye, Phys. Rev. D 24 ( 1994 ) 132; A. Martin, Phys. 
Lett. B 93 ( 1980 ) 338; C. Quigg and J. L. Rosner, Phys. Lett. B 71 ( 1977 ) 153; 
E. Eicheten et al., Phys. Rev. D 17 ( 1978 ) 3090. 

23. W. Kwong and J. L. Rosner, Phys. Rev. D 44 ( 1991 ) 212. 

24. Particle Data Group, L. Montanet et al., Phys. Rev. D 50 ( 1994 ) 1173. 

25. CLEO Collaboration, M. S. Alam et al, Phys. Rev. D 50 ( 1994 ) 43. 

26. M. Neubert, CERN-TH/96-55 (|hep-ph/9604412| ). 



27. C. Quigg, FERMILAB-Conf-93/265-T; C. H. Chang and Y. Q. Chen, Phys. Rev. 
D49 ( 1994 ) 3399; P. Colangelo et al., Z. Phys. C 57 ( 1993 ) 43. 

28. M. Bauer et al., Z. Phys. C29 ( 1985 ) 637; B. Grinstein et al., Phys. Rev. D39 ( 
1987 ) 799. 

29. A. J. Buras et al, Nucl. Phys. B 370 ( 1992 ) 69; Addendum, ibid, B 375 ( 1992 
) 501, B 400 ( 1993 ) 37 and B 400 ( 1993 ) 75; M. Ciuchini et al., Phys. Lett. B 
301 ( 1993 ) 263, Nucl. Phys. B 415 ( 1994 ) 403. 

30. C. S. Gao, J. L. Hu, C. D. Lii and Z. M. Qiu, Phys. Rev. D52 (1995)3978. 



12 



Table 1: Form factors in the SM and TCMLSM. f pen 9 and f anni represent form factors 
through electromagnetic penguin process and through weak annihilation process respec- 
tively. 



fi 


SM 


TCMLSM 


r-peng 
Jl 


-3.05 x 10" 10 


(-3.02 ~ -3.05) x 10- 10 


r-peng 
J2 


-1.52 x KT 10 


-1.52 x 10- 10 


fanni 
Jl 


7.10 x 10- 10 


7.10 x 10- 10 


fanni 
J2 


-1.70 x 10~ 10 


-1.70 x 10~ 10 



Table 2: The decay rates in the SM and TCMLSM. The T pen9 , T anni and Y total represent 
r ( B c — > D*7 ) through electromagnetic penguin process, through weak annihilation 
process and penguin + annihilation respectively. 



T(B C -> D* sl ) 


SM 


TCMLSM 




3.18 x 10" 19 


(3.12 ~ 3.17) x 10~ 19 


F anni (GeV) 


1.06 x 10" 18 


1.06 x 10~ 18 


T total (GeV) 


4.03 x 10" 18 


4.03 x 10~ 18 



Table 3: The branching ratio ( B c -> D*j ). The BR^f al and BRf gf LSM represent t h e 
branching ratio ( B c -> D*-f ) in the SM and TCMLSM respectively. 



r B c 


0.4ps 


l.Ops 


1.35ps 


UUSM 
nrL total 


2.44 x 10~ 6 


6.12 x 10~ 6 


8.27 x 10~ 6 


d tdTC'MLSM 
nn total 


2.44 x 10~ 6 


6.12 x 10~ 6 


8.27 x 10~ 6 



Figure captions 

Fig.l: The Feynman diagrams which contribute to the rare radiative decay B c — > 
-D*7 through electromagnetic penguin process. The blob represents the electromagnetic 
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penguin operators contributing to b — > 57, x 2 p and Xip are momenta of b and c quarks in 
the B c meson respectively, y 2 q and y\q are momenta of s and c quarks in the D* meson, 
respectively. 

Fig. 2: The Feynman diagrams which contribute to the rare radiative decay B c — > D*^ 
through weak annihilation process. In the SM, there is only W ± annihilation, in the 
TCMLSM, there are both W ± and vr^ annihilations. 

Fig. 3: VMD processes which contribute to B c — > D*^y with the vector-meson interme- 
diate states J/tp(^'). 
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Fig.2 
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Fig.3 
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